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Abstract 

We investigate electron transmission through a short wire with electron- 
electron interactions which is coupled to noninteracting metallic leads by tun- 
neling matrix elements. We identify two temperature regimes (a) Txondo < 
T < T'^^^^ = hvp/kBcl {d is the length of the interacting wire) and (b) 
T < Txondo T"^^^^. In the first regime the effective (renormalized) electron- 
electron interaction is smaller than the tunneling matrix element. In this sit- 
uation the single particle spectrum of the wire is characterized by a multilevel 
"quantum dot" system with magnetic quantum number S = Q which is higher 
in energy than the SU{2) spin doublet S = ±1/2. In this regime the single 
particle energy is controlled by the length of the wire and the backward spin 
dependent interaction. The value of the conductance is dominated by the 
transmitting electrons which have an opposite spin polarization to the elec- 
trons in the short wire. Since the electrons in the short wire have equal prob- 
ability for spin up and spin down we find G = G| + e'^/h < G < 2e^//i. In 
the second regime, when T — > the effective (renormalized) electron-electron 
interaction is larger than the tunneling matrix element. This case is equiva- 
lent to a Kondo problem. We find for T < TkoucLo the conductance is given 
by G = 2eV/i- These resuhs are in agreement with recent experiments where 
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for TKondo <T < T^'''^ the conductance G obeys e^/h<G < 2e'^/h, and for 
T < Txondoj G = 2e^/h. In both regimes the current is not spin polarized and 
the SU (2) symmetry is not broken. Our model represents a good description 
of the experimental situation for an interacting wire with varying confining 
potential in the transverse direction. 

PACS numbers: 73.21.Hb, 71.10.Pm, 72.25.Dc 
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I. INTRODUCTION 



Recent experiments in quantum wires show that the spin degrees of freedom combined 
with finite size effects give rise to new interesting effects in the baUistic transport regime. As 
an example we point out the new quasi-plateau with conductance G ~ 2(0.7 e^/h) observed 
by the Cavendish group [1]. A number of possible explanations have been suggested. Some 
of these explanations introduce spin polarization [2-4] and others are based on the Kondo 
effects [5,6]. Recently a number of groups have reported similar results. In particular we 
mention the results obtained by Reilly et al. [7,8] which have reported that the conductance 
varies in the range 0.5 — 0.7x2e^//iasa function of the electron density, length of the wire 
and temperature. 

In order to clarify this problem we investigate electron transmission through a short 
wire of length d with electron-electron (e-e) interaction coupled to non-interacting leads. 
Inspired by the experiments in Ref. [8] we consider a short wire of length d ~ 0.5-1 /im, 
which corresponds to a temperature T"'*^^ ~ hvp/kBd, where vp is the Fermi velocity and 
fcfi is the Boltzmann constant. We expect that the conductance will change drastically when 
the temperature changes from T < T"'*''^ to T < T"'*''^ [9]. 

We argue that the model of a short interacting wire coupled to non-interacting leads 
represents a good approximation to the experimental situation. In the experimental case 
the wire widens smoothly from the contact region to the reservoirs. In the contact region 
the wire is narrow, tuning the gate voltage allows a situation with only one propagating 
channel. In the reservoirs the wire is wide allowing for several propagating channels. (Due 
to the continuity condition from the reservoirs region to the contact region only one channel 
in the reservoirs is a pure propagating channel and the rest have a complex wave vector in 
the wire.) The presence of the e-e interaction can be separated into two parts, intrachannel 
and interchannel. Projecting out the channels with complex wave number we obtain an 
effective one dimensional channel. 

In the reservoirs the interchannel interaction renormalizes strongly the propagating in- 
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teracting channel. On the other hand, in the wire region the interchannel renormahzation 
is neghgible. As a result the reservoir region is described by a Fermi liquid contrary to 
the unrenormalized interaction in the wire region which is described by a one dimensional 
Luttinger liquid. The transmission coefficient between the two regions is described by the 
tunneling matrix element A which couples the one dimensional Fermi liquid in the leads to 
the one dimensional Luttinger liquid in the wire. The matrix element A is further reduced 
by projecting out the non-propagating channels. 

To solve this transport problem we construct the effective spectrum of the short wire. We 
find that the spectrum is built from charge-spin collective excitations and fermionic particle 
excitations. Using the renormahzation group (RG) we find that the spectrum of the short 
wire is controlled by the zero mode fermionic single particle states. The bosonic degrees of 
freedom are integrated out giving rise to an effective model for the short wire. As a result 
we obtain that transmission across the short wire is equivalent to transmission through a 
"multilevel state" at temperatures T < T""''^. The effective model describes the physics of 
length scale longer than d and contains new exchange terms generated by the interaction 
of the bosonic degrees of freedom of the short wire. We find that, since the conductance 
is given as the sum of the conductances for the different levels [see eq. (15b)], the S — Q 
state is higher in energy than the doublet S = ±1/2, i.e. the SU{2) doublet dominates the 
conductance and G ~ e^/ h. 

For temperature T — > so that T < Txondo < T""''^ the problem can be mapped into 
a Kendo problem which gives rise at T ^ to perfect transmission with a conductance 
G — 2e^/ h. This result is obtained only if the effective electron-electron interaction is larger 
than the effective tunneling matrix element. Therefore, this result is sensitive to temperature 
and length of the wire d. 

The plan of this paper is as follows. In the next section we present the model relevant to a 
short interacting wire coupled to two noninteracting leads. In Sec. Ill we obtain an effective 
model at (a dimensionless logarithmic) length scale I > Id = log{d/a), where d is the length 
of the interacting region and a is the inter-electron distance. In Sec. IV we present the 
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non- universal conductance at intermediate temperature and find that G ~ /h. Section V 
presents the computation of the conductance in the low temperature regime, T < Txondo, for 
which we find perfect transmission with the conductance G = 2e^/h. Finally, we summarize 
are main findings in Sec. VI. Calculational details are relegated to an Appendix. 

II. MODEL 

We consider a model of two non-interacting metallic leads coupled to a short interacting 
wire by tunneling matrix element, 

^ jjleads ^ jjwire ^ jj^ 

The method used in this paper is as follows: (a) The two reservoirs are described by two 
one dimensional non-interacting chiral fermions. (b) The one dimensional interacting chiral 
fermion describes a short one dimensional interacting region restricted to d ~ 0.5 - 1 /^m. (c) 
The transmission between the two regions is described by a matrix element A which couples 
the leads to the wire. 

Such a model represents the experimental situation in which an electronic waveguide is 
confined to the region —L/2 < x < L/2 and \y\ < D{x), where D{±L/2) = 1^ is the width 
in the reservoir and D{\x\ < d/2) = Dq is the width in the wire region with the conditions 
L ^ d and W > Dq. In the absence of e-e interaction this problem is solved within the 
Born-Oppenheimer approximation (see Ref. [10]). One finds that in the y direction we 
have a square well with transversal energies En{x) = {h^ /2m)[m^ /2D{x)Y' , n — 1,2,..., 
where n corresponds to the index of the channel. In addition, we have a negligible matrix 
element Zn,m{D{x)) which couples the channels. The value Dq determines the number of 
conducting channels Umax — 2kFDo/7r, where kp is the Fermi wave vector. For simplicity 
we consider the case with rimax = 1- This is obtained by tuning the gate voltage /x"'*''^. 
Due to the fact that the width D{x) varies in the x direction, the transversal energy acts 
as a one-body potential, En{x), and gives rise to backscattering. In the presence of a 
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ID Luttinger liquid we expect that the scattering potential En{x) will give rise to a low 
transmission coefficient between the two regions. [In the presence of e-e interaction the 
non-propagating interaction channels and the matrix element Zn,m{x) renormalize strongly 
the properties of the propagating channels in the reservoir. We know that for the case 
of coupled Luttinger chains the presence of small matrix elements between the channels 
is enough to give rise to a Fermi liquid. Therefore we expect that such a renormalization 
will take place in the reservoirs. Consequently the reservoirs can be described by a one 
dimensional Fermi liquid. In the wire region this renormalization is absent and we are 
left with a one dimensional interacting Luttinger liquid. The interaction will reduce the 
transmission coefficients between the two regions (Fermi liquid- Luttinger liquid).] Therefore 
we can describe the scattering matrix between the two regions by a matrix element A < 1, 
see Eqs. (6a) and (6b). 

The solution of our model, with leads coupled to a short wire, is obtained as a linear 
combination of the basis chiral operators CR^a, c^r^o- (right leads), CL,a, ^L,a (^^ft leads) and 
Xai Xa (short wire). The fermionic operators CR^a, CR,a'i CL,a, c'L,a Xa, xt constructed 
in each region separately for the case A = 0. The chiral fermion operators are obtained as a 
product of the zero mode fermion field and an exponential of the particle-hole bosonic fields. 
The bosonic field is periodic in each region. The boundary consitions are controlled in each 
region by the chemical potentials //J (left leads), /Ir (right leads) and /j,'^^'^^ (short wire). 

We construct an effective model for (a dimensionless logarithmic) length scale I > 
\og{d/a), where d is length of the interacting region d ^ 1 /xm and a is the inter-electron 
distance in the wire. Typical values of d/a are in the range 10^-10^. In the leads the inter- 
electron distance is different than in the wire but, since the leads are noninter acting, the 
Hamiltonian is scale invariant. Therefore, one can use the same inter-electron distance as 
in the wire. The only difference between the two inter-electron distances is incorporated by 
rescaling the coupling constant A. 

This construction is performed by using the Renormalization Group. At this length 
scale the bosonic degrees of freedom of the wire have been completely integrated out and 
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the fermion fields Xt{x), xti^) [^^- (2c)] have been replaced by the fermionic zero mode V^, 
[Eq. (2d)]. The short wire contains renormalized interactions gs, renormalized matrix 
elements and induced exchange interaction [e.g. Eq. (A6)]. The renormalized ratio X/gs 
and ^"'^''^ [Eq. (7c)] control the boundary condition at the interface between the reservoirs 
and the leads. For X/ gg < 1 we reproduce the Kondo boundary conditions. In the opposite 
limit we obtain the non-universal conductance G ~ e^/h. 

A. The leads 

The H^^°^' = H^ + H^ represents the left and right leads restricted to -L/2 <x< -d/2 
and d/2 < x < L/2, respectively. We will replace in each lead the fermion operators with 
chiral right movers, cl,c{x), c^L^ai^) (^^^ lead) and CR^crix), c^,o-(^) ("gl^t lead) [11]. This 
will be achieved by using open boundary conditions (OBC) for each lead. As a result the 
wave function of the electron vanishes aX, x = ±d/2 (i.e., the interface between the leads 
and the short wire) . The transmission between the leads and the short wire is described by 
a matrix element introduced at the interface between the two regions and is controlled by 
the chemical potential /j,^^^^. At the interface we will assume that the two wavefunctions 
in the two regions have an overlapping region of the order of the lattice constants. The 
transmission coefficient between the leads and the wire can be adjusted by varying the 
matrix element, the overlapping region and the chemical potential n^'^^^. (The concept is 
similar to the tight-binding model where localized states are used as a basis. By varying the 
hopping matrix elements we can describe the physics of free delocalized electrons.) 

We start with the chiral fermion representation of the electron operators in the leads. 
The electron field operators take the form: 

clA^) = e^'=^(^+'^/2)^^(^ ^ ^/2) - e-'^^^''^'^'^^ L^{-{x + d/2))- -L/2<x< -d/2. (16) 
Here cl,(j{x) is the fermion in the left lead expressed in terms of the chiral right moving 
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fermion L^{x). Similarly, for the fermion in the right lead we have: 

CR,a{^) = e'^^^'^-'^l^'^ R^{x - d/2) - e-^'=^(^-'^/^)i?^(-(x - d/2)); d/2 < x < L/2, (Ic) 

where CR^a{x) is the fermion in the right lead with the chiral right moving fermion Ra{x). 

Next, we perform the continuum limit and construct the Hamiltonian for the leads. We 
start with the Hamiltonian for the leads on the lattice (for each lead we have N sites): 



pleads ^ ^ [cl,,(-7Va-d/2)cL,,(-(iV-l)a-d/2)+...+4,(-a-d/2)cL,,(-d/2-£)+i/.c] 

-^0 E l4,,(iVa + d/2)c«,,((A^-l)a + d/2) + ... + 4,(a + d/2)c«,,(d/2 + £) + i/x]. {Id) 

Here a is the lattice spacing (~ inter-electron distance), e < a is the overlap region between 
the leads and the wire, is the hopping matrix element and Na — L/2 is the length of the 
lead. Equation (Id) represents the tight- binding Hamiltonian for the left lead with fermion 
operators c^^ o-, c|^^ and the right lead with c/^ o-, (^Ra- 

The OBC causes the boundary terms CL^a{—d/2) and CR^a{d/2) to vanish. This will occur 
when we substitute £ = in Eq. (Id). As a result we obtain h!^"-<^^ — jjieads^^^ _ p^j. 
£ 7^ we will have a transmission term from the leads to the short wire and in addition a 
non-zero boundary term. We express Hj^"''^^ in terms of the chiral fermions. 



H'^^^'^hvF Yl / dx[Ri{x){-id,)R,{x) + Ll{x){-id,)L,{x)]; hvp ^ 2toasm{kFa). 

(le) 

We also express Hamiltonian (Id) in terms of the even and odd chiral fermions 

IpeA^) = ;^[^'^(^) + '4'o,a{x) = -^[Ra{x) - L„{x)]. {If) 

Following the derivation in Appendix B we replace the leads Hamiltonian (Id) by 



pleads ^ pleads ^ ^^^^ ^^^^ 



/L/Z 
dx[^l„{x){-id^)^eAx)+^lA^){-idx)^o,a{x)]'^ = 2toasm{kFa). 

. -L/2 

{Ih) 



Hbc is the boundary term between the leads and the short interacting wire. This term 
takes a simphfied form in the Kondo regime where the "symmetric state" ipo,a{.x), ip\^{x) 
is screened out [the designation of the name symmetric or antisymmetric is given according 
to the original fermions CR^a{x) and cl,c{x)\- Hbc is replaced by the antisymmetric state 
'4^e,a{x), ipl^{x). In this limit we obtain for Hbc, 

Hbc ^ -2to{2sm{kFa)y ^ ^l,,(0)^e,a(0). (U) 

In this paper the zero-mode method [12] is used to solve the transport problem. Since 
the fermions in the leads are non-interacting, we will use the fermionic representation given 
by Eqs. (la)-(li). For the short interacting wire we will use the zero-mode Bosonization 
[see Eqs. (2a), (2b)]. For the sake of completeness we will also present the Bosonization for 
the leads. Because the leads obey L ^ oo the common belief is that the zero mode does not 
play any role. This belief is incorrect once we add particles to the ground state, in particular 
when we couple the leads to the external reservoirs [see Eq. (7b)]. For this situation it was 
shown in Ref. [13] that in order to study transport with Bosonization a nonvanishing DC 
current implies the presence of gapless modes, in analogy with the Goldstone theorem. The 
calculation in Ref. [13] is based on anomalous commutators. Similar results are obtained 
using the zero mode Bosonization [12]: 

For Lcr{x) we replace an^^ ctL,a, Nl,^, 6R^a{x) OL,a{x) and Xr^^Ix) XL,a{x), 

where au^^, OiL,a a-re the zero mode coordinates; Nr^^, Nl^^ are the number operators; and 
dR,aix), (^l,<t{x) are the particle- hole excitations. We have the following commutation rules, 
[o-R^ai Nr^iji] = iS„y, [q;l,o-, A^L,o-'] = The zero mode contribution to the Bosonized 

Hamiltonian is negligible since in the limit L oo, Nr^^/L and N^^^/L 0. 



9 



B. The short interacting wire 



Next, we consider the short interacting wire of length d<^ L. We will use the same Fermi 
momentum for the short wire as for the leads. According to the experimental situation we 
expect that the Fermi momentum in the wire region is smaller than the one in the leads (due 
to the change of the width) . This shift of the Fermi momentum in the wire is considered by 
adding to the wire Hamiltonian the chemical potential ^"'^''^ [see eq. (7c)] which shifts the 
wire ground state. 

The finite size of the system and the boundary conditions for the short wire allow us to 
introduce zero mode excitations. The fermion da{x), dl{x) on the short wire is represented 
by anti-periodic right chiral fermions Xct{x), 

d,{x) - e''-'xa{x) + e-''-'xa{-x). (2a) 

Here da^{x = d/2) = d„{x = —d/2) = is obtained by demanding anti-periodic boundary 
conditions for the chiral fermion in the short wire. 

Xa{x + d)^Xa{x)e'\ (26) 

Using the zero mode Bosonization introduced in Ref. [12], we obtain the representation 

Xa{x) = i/^e(^''^/'^)(*''^-V')^e^^'^'^(^\ (2c) 

= ^e-^'^- = <T=T,i- (2d) 

Here andp^ are the zero mode coordinate and "momentum", \(ia^'Pa'\ = i^aa'^ respectively. 
The "momentum" — 0, 1,2, ... measures the change of fermion number with respect to 
the filled Fermi sea. (\)cr{x) is the non-zero mode (the particle-hole) bosonic excitation. The 
operators yj, Va act as creation and annihilation fermion operators, Vj|pcr) = \Va + !)• The 
operators V^J, V„ obey anticommutation relations, {yj, = ^^aa'i {Yci ^J'} — Kr'} — 
and commutation relations with the number operator p^.; [Po-; Kr'] = ^<t<t'^ct'; b^; Kt'] = 
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The form and condition in Eqs. (2a)- (2c) are chosen for the following reasons: (a) We 
intend to construct the bosonic representation of the chiral fermion field xAx) which is 
given as a product of the zero mode fermion exp (i27r/(i)(po- — 5)x and the exponential of 
a periodic bosonic field (t)a{x). The boundary condition of the fermion field da{x) in Eq. (2b) 
is determined by the value of 5. In Eqs. (2b) and (2c) we have 5 — 1/2. The representation 
in Eq. (2c) is given for 5 = 1/2 and has no fermions in the ground state, (po-) = 0, Po- =: Po- 
where : Po- : stands for normal order and (po-) nieans the expectation value with respect to 
the filled Fermi sea. 

Equation (2c) is obtained in the absence of a chemical potential of the wire. Once a 
gate voltage 2Vg is apphed, we must add to the Hamiltonian a term //"''''^(p-f + pj [see 
Eq. (7c)]. The presence of the nonzero chemical potential ^^'''^ shifts the ground state 
occupation. Consequently (po-) ^ and p^. — 1/2 is replaced by : p^- : +(P(t) ~ 1/2. The 
shift (po-) — 1/2 can give rise to a change in the boundary conditions. We observe that 
the coupling Hamiltonian in Eqs. (9d) and (11c) below is controlled by Po- = (P(t)+ '■ Va 
with (po-) determined by the renormalized chemical potential E„{d) [see Eq. (9h)]. In the 
intermediate range of temperatures (see Sec. IV) the effect of the chemical potential in the 
wire is considered by taking the expectation value with respect to p^- [see Eq. (15b)]. 

The physical reason for the representation in Eqs. (2a-2c) has to do with the fact that 
at the (dimensionless) length scale I > log{d/a) the short wire is replaced by an impurity 
atom which in the limit T — > corresponds to a single impurity. In this regime A <^ 1 and 
the e-e interaction Qg is large (see Sec. V). This gives rise to the Kondo physics. 

Using the chiral representation given by Eq. (le), we obtain the Hamiltonian for the 
short wire Hubbard model 

jjwire jjwire _j_ jjwire _j_ jjwire (3o) 

where 

rd/2 
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(36) 



and U is the original Hubbard interaction, U = 2U and vp = 2vf- The Umklapp and the 
backward term Hamiltonians are respectively given by [14]: 



1 _ fd/2 , , 

nr^^-U dx{x\{x)x\{x)xi{-x)xA-x)e''"'-^ + H.c.}, 

Z J—d/2 



(3c) 



1 ~. f'^/'^ 

l-d/2 

-U dx{x]{-x)x\{x)xi{-x)x]{x) + H.c}. {3d) 

2 J-d/2 

Next we bosonize the Hamiltonian given in Eqs. (3a)- (3d). Hq'^'^^ can be replaced by 
ifi"=°) (the zero mode) and Hq""^^^ (non-zero mode part): 



TTWire 



TTWire _ rrin^G) . Tjin^G) 



(4a) 



Here H^f^^ represents the charge density and H^f^^ is the spin density. 



HiT^ - / ' dxVs{dM^))\ 

' -d/2 



d/2 
d/2 



(46) 



(4c) 



(Ad) 



The renormalized bosonic fields, (j)c{x) and (f>s{x) (charge, spin) are related to (pdx) and 
(j)s(x) by: 



K 



-1/2 



K 



1/2 



(t^cis){x) = (0cw(a;) + (Pc{s){-x)) + — ^ (0c(.)(a;) - 0c(.)(-a;)) , (4e) 



with Kc and i^^ given by: 



1-U/71V+ , U 



\ l + C//7rv+' 



l + U /7rv_ 



= {ii^ ( 1 — 



?7 



and charge and spin velocities 



Vr = 
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The zero mode part of Hq allows viewing the short interacting wire as a multilevel 
"quantum dot" with single particle energy given by: 

where p^- — 0, 1, 2, ... represents the additional charges with respect to the filled Fermi sea 
\p-l = 0,pi = 0). Typical values considered in Ref. [8] are do ~ 0.5 /im. This corresponds 
to energies eo = hvp/^d — 10"^^ (vf/c) {do /d) Joule, since vp/c ~ 10~^, eo is of the order 
of meV or temperature T'"'*''^ — eo/ks — 7.6 x 10^{vF/c){do/d) which corresponds to a few 
Kelvins. 

The Bosonization of the Umklapp and backward terms involves both the bosonic degrees 
00- (x) and the zero modes p^j, and fermionic operator V^. 

27r 



^wire ^ J^^^^ i^__v\v^vtv^ cos [(4fc^ - G)x + ^(pt + Pi) + Vs^iM^) - M-^))] ] 

(56) 

We find that the Umklapp term is controlled by the charged boson (f)c{x), the charge of the 
"dot" Pc — Pi +Pi and the reciprocal lattice vector G — 27r/a. This term is highly sensitive 
to the electron density, namely kp = {n/2a){l — 6). For density S and wires of length d 
which satisfy {Akp — G)d> 27r or 6{d/a) » 1, the Umklapp term can be neglected. Since 
d ~ 0.5 - 1 /xm the deviation from half filling must be of the order less than 1%. Under this 
condition the Umklapp term will not give rise to a charge gap. 

Next, we consider the backward term, Hq"^^^. In the limit d — > oo this term renormalizes 
to zero driving Kg ^ 1. For finite d this is not the case: 

27r 



TTWire 



J-d/2 "^"^ I ^^^i [J^^^ ~ ^^^^ + Vs^iM^) - M-^))\ I ■ (5c) 



The background term is controlled by the bosonic spin density 0s (x) and spin excitations 

Ps=Pi - Pi- 
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C. The coupling Hamiltonian 



The coupling Hamiltonian is given by 



Ht^ Yl HLcl,i-d/2 - e)d,{-d/2 + e) - tRC^R^,{d/2 + e)d,{d/2 - + H.c]. (6a) 



In Eq. (6a) we have introduced an "overlap" £ ~ a between the electrons in the wire and 
leads in order to allow transmission. (Due to the boundary condition when £ = we have 

Ht = 0). Next we assume ti = = t and make use of the chiral representation [see Eqs. 
(lb,lc,lf,2a)]. We find: 

a=1,lR=±(,d/2-e) '^/'^ 

The tunneling matrix element A = A(A)°, where A = 1/a is the cutoff and A = 
2y/2tsm{kFa) = {t/to)EFSm{kFa), Ep = Uvp/^a. 

The matrix element A obeys A < 1 (this being a result of the e-e interactions in the leads 
which have been integrated out). This value of A is fixed by the transmission strength t/to 
and the overlapping region of the wavef unction £ ~ a. We considered the short wire as a 
localized state which is coupled to the leads. ControUing the strength of A we can describe 
extended solutions in spite of the fact that we have started from a localized picture (this is 
the philosophy of the tight-binding method) . 

We argue that for our problem this is a good starting point since renormalization effects 
further decrease A and give rise, for the (dimensionless) length scale / > log((i/a), to an 
impurity problem which can be mapped into a Kondo problem for T ^ 0. 



D. Computation of the current 



To compute the transmission current we add the reservoir Hamiltonian Y: 




(7a) 



leads 



(76) 
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with (1/2) (//^ + — (1/2) (/x^ + /x^) = eVos the voltage difference, and Nl^o- — 

L /2 + L l'2. + 

l_j^l2dxL\{x)L„{x)^ Nr^cj = j_i^i2dxR\{x)Rcr{x) are the fermion densities in the leads. 

^we^^we^^^^^^)^ (7C) 

with jjL^'^'^'^ — eVa being the gate voltage apphed to the wire. The term y^'^^ in Eq. (7c) 
allows for a nonzero number of fermions with respect to the leads Fermi sea. When /x'^*''^ ^ 0, 
kp in Eqs. (2a)-(2c) is shifted to kf^'^'^ = kp + {27r/d){pcr) and p^r is replaced by : p^. :. This 
procedure accommodates the experimental situation, kp'^'^^ < kp (//"'^''^ can be tuned to 
(Pa) < needed in the experiments). We note that the conductance formula given below in 
Eq. (15) is sensitive to this tuning. The current operator /o- is given by 

^'^ R=±{d/2-e) 

The non-equilibrium value of the current is obtained by performing the thermal expec- 
tation with respect to the grand canonical Hamiltonian, H — H — Y , where H is given in 
Eq. (la) and Y in Eq. (7a). As a result we obtain 

Equation (7e) will be used below for the computation of the current. 

Equation (7e) with the reservoir Hamiltonian Y^^°^^ = eVos^hl'^ - ^VdsNr/^ [see Eq. 
(7b)], with iiL = eVDs/2, iiR = -eVDs/2 and Nl = + N^i, Nr = Nr^ + A^^^ allows 
to compute the current taking into account the boundary imposed by the reservoirs. Con- 
sequently, the density in the right and left leads obey the boundary conditions given in Eq. 
(14c). This method has been used in Ref. [13]. In particular, it has been shown in Ref. [13] 
that a one-dimensional interacting fermion (with no backward and Umklapp interaction) 
coupled to non- interacting leads has a universal conductance G — 2e^/h. This result has 
been shown to follow from the dynamical requirements, [Nr, Ni\ = [H, Nr] = [H, N^] = 0. 
{H is the interacting Hamiltonian constructed in Refs. [15-17]). The non-renormalization of 
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the conductance by electron-electron interactions has been shown originally [15-17] to be a 
consequence of the strong influence of the boundary conditions imposed by reservoirs. The 
same results have been obtained in Refs. [12,13] using the y'^'^'^* Hamiltonian [see Eq. (7b)]. 
The advantage of the formalism used in Ref. [13] and the zero mode calculation given in 
Ref. [12] is the simplicity. On the other hand the result derived in Ref. [15-17] are based on 
the exact integration of the electrostatic potential V(x) along the one dimensional channel 
/-L/2 where p{x) is the electronic density. For problems with backscattering 

[induced by the one body potential En{x)] and backward interaction, the method given in 
Ref. [13] is easy to use. The simplicity of the reservoir Hamiltonian Y'-^'"^'^ = urNr + iilNl 
(see Ref. [13]) causes only the zero mode coordinates (or the Goldstone modes) to be af- 
fected by the external boundary conditions. At flnite temperature we use Eq. (7e) and at 
T = 0, Y^^"''^^ is part of the Hamiltonian which will affect the zero mode variables. When 
yieads _^ Q ^^le grouud statc is shifted, consequently the zero mode coordinates become 
o;r,<t(^) (^R,ai't) ~ eVost/'^fi and aL,o-(^) <^L,ai't) + eVDst/2h. As a result the chiral 
fermion in the leads will be modified. The current will be computed by replacing the tun- 
neling Hamiltonian given in Eq. (6b) by a time-dependent Hamiltonian. Therefore, Ra{x) 
and L^(x) in Eq. (6b) and (7d) will be modified: R^{x,t) — > exp(ieVD5t/2^)i?o-(3^, ^) and 
La{x,t) exp{-ieVDst/2h)L^{x,t). 

The boundary condition affects only the zero modes OiR^a, oll^^ (leads) and (wire). 
Therefore, when we integrate short distance behavior given by (t)a{x) in Eq. (2c), this has 
no effect on the zero modes and thus on the boundary condition. According to Refs. [15-17] 
this does not appear to be correct since the result depends on the integral of the electrostatic 
potential V(x) along the one-dimensional channel. The resolution of this problem can be 
achieved by replacing the electrostatic potential V(x) with an averaged electrostatic poten- 
tial. Using this approximation will allow us to neglect the effect of the boundary on the 
short distance renormalization processes. For the rest of the paper we will work at finite 
temperature. We will use Y^^"''^^ — iirNr + UlNl according to Eq. (7b) and will compute 
the current using Eq. (7e) in agreement with Ref. [13]. By doing so we ignore the short 
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length fluctuations of the electrostatic potential considered in Refs. [15-17]. 

III. EFFECTIVE MODEL 

The effective model at the (dimensionless) length scale, I >ld = ^og{d/a), T < hvp/kBd 
can be mapped to an impurity model. To obtain the effective model we integrate out the 
short wire degrees of freedom. The wire field operators are given by 

This operator contains bosonic degrees of freedom (/>(r(x), which can be integrated out as 
regular bosonic fields, and a discrete fermion number which measures the added changes to 
the wire, — 0, 1, 2, ... and V^, are the fermionic operators (see Eq. (2c)). 

This integration is performed in two stages: (a) Integration of the bosonic degrees (t)a{x). 
This step is performed with the help of the Renormalization Group. As a result of this KG 
the short wire Hamiltonian and the tunneling matrix element A are renormalized according 
to the sine-Gordon scaling equations. In addition, new magnetic exchange terms are induced 
by the RG. After performing the RG we use experimental considerations, namely that the 
studies of conductance have been performed below a few Kelvins for wires with length of 
the order of 1 ixm. This means that the relevant physics occurs at a (dimensionless) length 
scale I > Id = \og{d/a). Consequently we compute the relevant coupling constants at the 
scale I — Id- We substitute A = 1/d as our new cutoff and rewrite the model in terms of the 
scaled coupling constant and new fields, 

Xc{R = ±((^/2 - e)) ^ v;e'(2Vd){P.-i/2)ii_ (g) 

Note that in Eq. (8) the bosonic field has completely disappeared. 

The model which is obtained for the short wire corresponds to a multilevel "quantum 
dot" system characterized by the quantum number p^. Equation (la) is replaced by: 

H{d) = H^^'^'^'id) + H'^'^^d) + Hrid) + H\\{d) + H^{d) + Hp-p{d) + H^{d). (9a) 
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Here, W^"''^^{d) is the same as H^^"^^ except for the fact that V'e.o- — '4^e,a and ■00,0- — ^ V'o.o-, 
where ^J^, z/'e,^, ^^g. and t/'o,cr represent the fermions in the leads with the cutoff A — 1/d 
instead of A = 1/a. 

fjwtre^^-j represents the renormahzed short wire Hamiltonian and is given by: 

In Eq. (9b), rj — 2U {Any / deo ~ 1 is the dimensionless interaction parameter; gdld) and 
Qsild) are the renormahzed Umklapp and backward interactions, respectively, and p„ repre- 
sents the quantum number of the multilevel system. When {ikp — G)d > 2t: we can neglect 
the Umklapp term. On the other hand the backward interaction term vanishes in the limit 
d ^ oo. 

The value of gs{ld) is determined by the sine-Gordon scaling equations [14]. For a finite 
d and Kg > 1 we find. 

Since -^^5(0) > 1 and d is finite, Ks{ld) > 1 and therefore the critical scahng d — > 00, — > 1 
with gs{l 00) ~ ^s(0)/(l + fi's(O)/) is not applicable. 

We also approximate the Hubbard interaction t/ by a screened Coulomb interac- 
tion with a dielectric constant k, — s/sq ~ 1. We approximate U/a ■2:^ /{Ansi^jKa. 
As a result we find that U is related to the energy of the short wire eo = hvp/d, 
U ~ ^^{c/vF)d, where c ~ 3 x 10^ m/sec is the speed of light, c/vp — 10^ and 
K ~ 10. Using these values we estimate -^^^(O) = ^^(1-1- C//7rtof-)/(l — U/irtoV^) ~ 
^(1 + {l/137){c/vF){l/7rK))/{l - (l/137)(c/wir)(l/7rK)) ~ 1.5. This estimate allows us to 
replace Eq. (9c) by gs{ld) ^ {€o/2Ti^K){l/137){c/vF){d/aY^-^^^\ For d/a ~ 10^ we find 
that gs{ld) ~ ^0) which is a few meV. 

The next term is the coupling Hamiltonian Hrid), which replaces Eq. (6b) by 

/-L/2 ^ ^ / 27r \ 

Hrid) ^2iXY. / dx5(x)[iljl,{x)V,cos(kF{d/2 - e) + — (p, - l/2)(d/2 - s) 
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- COS (kF{d/2 - + ^(p, - l/2)(ci/2 - £)) V^JoA^)], {9d) 

where A = {2/y/7rd)tsm{kFa) — AA^/^, A = 1/d and A = (Vv^) sin(A;i;'a). This result was 
obtained within an RG calculation; see Appendix A. 

The next three terms represent the induced magnetic interaction and effective impurity 
energy: 

H\\{d) = 2 Y: + I\\{ld) cos (2^^(ci/2 -s) + ^(p, - l/2)(d/2 - e))] V.t ^(o)4,,(0)V;tV;. 



(9e) 

The coupling constants J\\ and /|| are given in Appendix A (see Eqs. (A7,A8)). The transver- 
sal exchange term is given by 

/27r A 271 



H^{d) = 2J^{U) 



cos (^-^Ps{d/2 - e)j + cos (^{2kF + -^Pc){d/2 - e] 



x[<T(0)^oan^VT +<^(0)4,T(0)V^/V^i]' (9/) 
where Ps — — Pi and Pc — P^ + Pi- The coupling constant J±{ld) is given in Appendix 
A (sec Eq. (A9)). Wc observe that the coupling Hamiltonian in Eq. (9d) and the induced 
interactions in Eqs. (9e)-(9g) are sensitive to the value of the overlapping region e and 
chemical potential /i"'*''^. The use of £ 7^ and ^""''^ ^ changes the boundary condition 
to S ^ 1/2 [see Eqs. (2a)-(2c)]: {p„ - l/2)(rf/2 - e) = {: p^ : -5){d/2 - e), therefore 
(5 = l/2-(p,). 

The induced two particle term is given by Hp-p{d), 

Hp-p{d) = Mk) cos [^Ps{d/2 - £)) + cos (^{2kF + '^Pc){d/2 - e)j 

[i^i^mlMViVi + y/y/4a(o)4,T(o)]. (9^) 

Hp_p{d) gives rise to two particle transmission at low temperature. The two particle term 
dominates when the single particle transmission controlled by A vanishes. 

Equations (9e) and (9f) represent the induced magnetic interaction. These equations are 
obtained from the first two terms given in Eq. (A6). The particle-particle term given in 
Eq. (9g) is obtained from the last two terms in Eq. (A6). The explicit form in Eq. (9e) is 
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obtained once we integrate out completely the bosonic degrees of freedom (f)a{x) of the short 
wire. At scale / > \og{d/a) we replace in Eq. (A6) fermion fields Xa, xl- by the zero mode 
fermions V^^V^. 

The last term in Eq. (9a) is the induced single particle energy 



where 



Eaid) 



1 y r J . c . 27r 



J||(/d) + mid) cos ( {2kF + -jPa){d/2 - e] 



>0, 



and ~ -E'o-(rf) are obtained after performing the normal ordering in Eq. (A6). When a 
gate voltage is apphed, Ea{d) is shifted by the chemical potential of the wire yu"'*''^ [this will 
be the case in Eq. (15f) below where the single particle energy ei is shifted, ei — > ei — eVc]- 
The set of Eqs. (9a)- (9h) complete step one of the renormahzation group. 

In the second step we map the problem to an impurity model. This step is performed 
using the impurity Hamiltonian (Eq. 9b). We have to project out high energy states: 
we project out the states with p^. = 2, —1, —2, .... As a result of this projection we keep 
the same structure as given in Eqs. (9a)- (9h) with the difference that p„ is restricted to 
Po- = 0,1,(7 =t,i and the coupling constants J\\{ld)-, -^||(^d) and J±{ld) are replaced by 
(J||(Zd)) = J||(Z<i), (7||(Zd)) = 7||(Zd), and {J±(ld)) = J±(ld) with "(...)" standing for the 
averages. Prom Appendix A we find that in the limit d ^ a the induced exchange coupling 
constants J\\{ld), I\\{l'd) and J±{ld) vanish. The projected Hamiltonian h replaces H{d), 
H{d) h 

h = H^^'^'id) + himp + /iT + h\\ +h^ + hp, (10) 
where H'^^'^id) + H^{d) himp, Hrid) hr, H\i(d) /i,,, H±{d) h±, and Hp^p hp. 
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IV. CONDUCTANCE IN THE INTERMEDIATE TEMPERATURE REGIME: 



For Txondo <S T" < hvp/kBd ~ J'""''^ we use the impurity model obtained in Eq. (9b). 
In particular, the impurity Hamiltonian with = 0,1 and a =t, i is given by, 

himp = + + WJ + gs{ld)V^^W;V^Sp^.p^,o - E,{V^^V^ + y/Vl), (11a) 

where gs{ld) is the renormalized backscattering term given in Eq. (9c). [We have neglected 
the Umklapp term (Akp — G)d > 27r]. In Eq. (9h) we replace Er^{d) by ~ {J\\)/'Kd. 

In this temperature regime we ignore the two particle transmission and replace h\\ + /i^ + 
hp-p by hf 

h\\ ^ (ill) E i^l,Mk,Mv^y^, (ii&) 

with (J||) ~ (2t^/T:vF)\og{d/a). 

The transmission term in Eq. (9d) is replaced by hr- 

hr ^2A ^ cos(7rp,)[<,(0)V; - Vji^oA^)]- (He) 
^=1,1 

We have approximated the cosine term in Eq. (9d) by cos(7r|Jo.)- This approximation has a 
negligible effect on the renormalized tunneling matrix element A. The matrix element A is 
fixed by using e 0, S — 1/2 or equivalently e — 0, S 1/2, see Eq. (2c). 

The impurity spectrum is characterized by the ^7(1) x SU{2) Kac-Moody primary states 
of the form \p],Pi). When pc = Pt is even the state |pi-,Pj,) = \pc = even,Ps = 0) is an 
SU (2) singlet with Ps — Pi — Pi — 0. On the other hand, when Pc — Pi + Pi is odd the state 
\PhPl) ~ \Pc — odd,Ps ± 1) is an SU{2) doublet Ps — ±1 and spin Ps/2 — ±1/2. 

The Heisenberg equations of motion are obtained from the Hamiltonian in Eq. (11a) 
with the zero mode anticommutation relations {V^, Va'} — 2Sa,a', V}'} — {Va, Va'} — 
and commutation relations [Po-,V^/] = ^ayVj', bo-jK-] = — (^<T,cr'Kr'- 

We replace the number operator p^- by V^V^ = 2Ti{d/a)p„. The factor 2Ti[a/d) is related 
to the different renormalization of p„ and V^V^. Next, we determine the eigenvalues for the 
impurity model [Eq. (11a)], 

himp\Pl,Pl) = £{Pl,Pl)\Pl,Pl)- 
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We find that due to the fact that gs{ld) > eo and Ed > 0, the SU{2) doublet states \pc — 
odd,Ps = ±1) have a lower energy than the singlet state \pc = even,Ps = 0). Therefore we 
expect that the current will be controlled by the SU (2) doublet. 

In the remaining part of this section we compute the current which is obtained from 
Eq. (7d) with Ht replaced by hr [see Eq. (11c)], 

la^Y cos(7rp.)[4V(0)V; + V;t4,.(0)]. (12a) 

To find the current from Eq. (11) we have to solve for the scattering states [18,19]. The 
scattering states will be obtained within the Heisenberg equations of motion for V^, il)e,a 
and ipo,a- 



r 1 ^ 1_ 
ihV^ = [V{, himp + hT + /i||] = -eo PT - 2 + + 2^sPi^p^-pi, 



+EaV; - i(2A) cos(7rpT)^^o,T(0) + {JM^momVv (12&) 
ihiP^-^ix) = [4t(0), himp + h + H^^'^^d)] = hvF{-id^)'iPoA{x) 

+J\\P^'4Jo,^{x)6{x) + i{2X) cos(7rp|)V/(5(a;). (12c) 

The zero mode p^ which measures the charge fluctuations with respect to the Fermi 
energy obeys the Heisenberg equation of motion 

ihpa = \pa, himp + hr + h\\] = \pa, hr] = -2Acos(7rp^)[V'i^(0)V; + V;^4,a(0)]. (12e) 

In Eq. (12b) we have used the notation: (1/2)^^, = gs{ld)27c{d/a)(l/€o) and J\\ = 
{J\\)2n{a/d). From Eq. (12c) we find that ipo,a{0) scales with 2A and the spectrum of 
Va is controlled by the single particle energy eo. Therefore, we conclude that we can neglect 
the time dependence of p^ when the single particle energy is larger than the transmission 
energy, (2A)^/eo < 1. In this limit we can neglect the fluctuations Spa{t), Pa{t) — Pa + ^Pait) 
and we make the approximation (pait)) — Pa- 

Next we use the bosonic representation of ipo^a and replace, ■?/'j,-|.(0)'?/'o^-|-(0) = (l/27ra) + 
{l/2'!r)dx0o,i- Here 9o,i is the bosonic field in the leads and is the renormalized impurity 
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energy, Ea^ Ed+ (J||)(Vit (0)4(0)) = + {J\\) /'^T^d. We neglect the dJo^iO) term in Eq. 
(12b) and solve Eqs. (12b) and (12d) together. We expand ipo,a{x,t) in scattering states 
Ue,<j{x) with the eigenvalues E. Then 

E 

K{t)^Y.M'rit)Va,E, (136) 

E 

with AE^ait) = AE,ae~^^^^^^^ ■ We find from Eqs. (12b), (12d) and the approximation 
{5pa{t)) — that V*^E is proportional to C/b,(t(0), 

y* ^ i2Acos(7rpt)?7j;,|(0) 

^'"^ 6obT-| + iPi + |pAr-n,o]-^<.-£^' 

and 

{E + hvF{id,))UE,^{x)^JllP^S{x)UE,^{x)+i2Xcos{np^)V^^^^^^ (13d) 
We substitute Eq. (13c) into Eq. (13d) and obtain the scattering state solution Ue,i{x): 



l^i-x) + i^ix) ' 



(14a) 



l + i{2Xy/Sj\ ' 

where //(x) is the step function ii{x > 0) — 1, ii{x < 0) — and S is the resonance energy 
function 

S = S{E-p^,p^) = eo(pT - 1/2 + (r;/2)pi +^,Pi5^^_p^,o) -E^-E, (146) 

where = hvp/^d, rj = ?7(47r)^/eo(i, and is given by Eq. (9c). According to Eq. (12d), 
'ipe,a is given in terms of plane waves, i^e,a{x) — YIiE BE,a{^/ 'sfL)e^^^^^^'"''^^^ , and ^o,a{x) — 
J2E^E,aUE,a{x), whcrc Ue,(t{x) is the scattering state (Eq. (14a)). 
In the leads the presence of the voltage difference Vds gives 



{{RIR.)e) = {a^E,aaE,a) ^ fpD [e + ^) , 

{{LlLM = {h'^EM^) ^ {e - ^) , (14c) 
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where fpoiE) is the thermal Fermi- Dirac function. The right and left moving operators 
ttE^a s-nd hE,c s-re related to ■0e,(T ^-nd -^o^c, by the relation: 



B 



1 



{aE,a + bE,a)- 



{Ud) 



We substitute into Eq. (12a) the Eqs. (13a) and (13b) using the eigenfunction UE,a{x)- 
We take the thermal expectation with respect to the leads reservoir [see Eq. (7b)] and find: 



(/t) = ^cos(7rpt)gi [(i?t._^yls,t)VT,i. + (4,T^i^,T)V'T: 



le 



{2\f 



Ah 



' ? S{E) 



h 



dE 



cos(J||Pt) + ^ I , :^/J sm(J||Pt) 



(2A)4 + 52 



X 



Jfd ( ) — fpD ( -E" + 



S \S^ + {2XY 
eVns 



2 J -v- ■ 2 /J 

In Eq. (14e) 5" is given by Eq. (14b). When we compute the current (/j^) we have to replace 
PT ~^ Pi Pi ~^ Pi -^1- (14e). Next, we perform the expectation value over the impurity 
states at temperature 1//3 using the result in Eq. (14b): 



((4)) 



(15a) 



Expanding the Fermi-Dirac function in Eq. (14e) around E = Q gives for the current the 
conductance: 

{{I^))_e'// (2A)^cos(J||p^) \\ 



e 



Vds h\\(2X)^ + (S{p^,p^)y//' 
where S{pi,pi) = S{E — 0;p-f,pj. Here ((...)) stands for the thermodynamic sum over 

the impurity states: \p^ = l,pi = 0), \p^ = 0,pi = 1), \p^ = 0,pi = 0), jp^ = = 1)- 

Equation (15b) is dominated by the resonance term S{pi,pi) which obeys (2A) > \S{p-\,pi)\. 

In obtaining Eq. (15b) we have neglected the term proportional to sin(J||p-f) in Eq. (14e). 

This approximation is justified given the fact that J\\ <S 1, see Eq. (A14). In order to 

perform the expectation value in Eq. (15b) we consider the even and odd state eigenvalues 

e{pi,Pl) and the resonance term S{pi,pi). 
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The even states: 

£(pt = o,Pi = 0) = 0, s{p^ = o,Pi = 0) = - (I + 

£(PT = = 1) = 2ei + C/e-e, = = 1) = S{p^ = 0,pi = 0) + 2ei + C/e-e- (15c) 

The odd states: 

£(PT = I'Pi = 0) = ei, = l,pi = 0) = ei; 

£(PT = 0>Pi = l) = ei, S{p^ = Q,p^ = l) = e, + (^-l]eo, (15d) 



.2 

where ei = (eo/2 — E^^) is the "impurity" single particle energy which is below the Fermi 
energy, t/g-e = £o(^/2 + gs) is the effective "Hubbard" interaction and S(j>^ = 0,pi — 
1) — = l,pi = 0) = (77/2 — l)eo is the effective spin polarization energy. The energies 
^(PhPl) control the thermal expectation values. We observe that due to the electron-electron 
interaction and the single particle energy the odd states, e{pc — odd,Ps = ^1), can 
have lower energy than the even states. 

The denominator function S{p^,pi) controls the quantum weights for the different states 
\pl,Pl) and it is asymmetric under the transformation p^ —>■ Pi, Pi —>■ Pi- For the odd 
states it has the property for 7-f, S{pi — 0,pi — 1) S{pi — l,pi — 0). [Similarly, for I^, 
S{pi — 0,pi = 1) ^ S{pi = l,pi = 0)]. Due to these properties we note that the current is 
dominated by one of the two states. We use Eqs. (15b)-(15d) to find: 

^ ^ e\ (2A)'- (2A)^-os(J||) (2A) ^-os(7||) 

^ h ^(2A)4 + (-5(0, 0))2 (2A)4 + {S{1, 1))2 (2A)4 + {S{1, 0))^ 

(2^)' -e-/^^(0'i)][l + e-^^(^'^) + e-^^(^'°) + e'^^^'^'^^J-^ (15e) 



(2Ar + (5(0,l))2 

The conductance in Eq. (15a) depends on the effective single "particle" energy, ei = 
eo/2 — Ed < 0, the vacuum energy, ec = eo/2 + E^ = —S{pi — 0,pi — 0), the effective two 
particle energy Uee = eo(^/2 + ^s) > and polarization energy tg = eo(^/2 — !)• In the 
presence of the gate voltage Vq the single particle energy is shifted ei ^ ei — eVc- (The 
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single particle energy ei is controlled by the renormalized wire energy and by the chemical 
potential shift fi^'''^ = eVc-) 

Next we introduce dimensionless parameters: (2A)^/eo = T, ei/eo = ei, eVc/eo = Uq, 
C/ee/eo = Uee, Cc/co = ^c, Cs/co = ^s, eQ/ksT = T'^'^'^^/T, and cos(J||) ~ 1. As a consequence 
we find that the conductance takes the form: 

(7, - ^ r I -^[2(ei-C7G)+(r,/2+flJ]. 



2 



/i Z ^2 + + jj^^2 r2 + [2(ei - Ug) + (77/2 + ^J] 

^^-m-ua) ( . + J! Y (15/) 

where Z = l+e-'^P{ei-t>G)+{r?/2+s,)]_^2e-'^('"i-^G). Due to the repulsive interactions 2ei+t/ee > 
and ii < 0, we find that for T^*''^ > T the even part Pc = 0, 2 in Eq. (15f ) can be neglected. 
Consequently for Ug — Eq. (15f) is replaced by: 

,2 1 



+ 



(15^) 



P + ef r2 + (ei + r;/2-l)2_ 

We investigate Eq. (15g) for ei < 0. As a result the term e^^'""^''/^)^! — > when 
rj^wire jrj^ ^ -Q^g repulsive interaction ri/2 we have a situation where one of the two 

terms satisfies the condition (ci/F)^ < 1 or (ei + 7y/2 — l)^/r^ < 1. Choosing ii+rj/2 — 1 — 
we find: 

1 + 



^ 2/i 



e2 



r2 + {rj/2 - 1)2 

2 2 

A similar calculation for gives Gj^ ~ Therefore G = G-f + Gj^ ~ ^. This result 
is not based on a spin polarized state but makes use of the asymmetry betweed odd and 
even states such that only one state gives a resonant contribution (2A)^/((2A)^ + S"^) ~ 1 
for ,5(0, 1) or ,5(1, 0) but not both. 

The result obtained above depends on the number of electrons in the wire. We consider 
the case Pi +Pi — Pc — 1- (The even case Pc — 0,2, ... can be ignored since these correspond 
to high energy states and therefore give neghgible contribution to the current). When Pc — I 
we have one electron in the highest state in the wire. Suppose that an additional electron 
transmits into the wire. Here we can have the following two situations: (a) The spin of the 

26 



electron in the wire is the same as the spin of the transmitting electron; (b) The spin of 
the electron in the wire is opposite to the one which is transmitting into the wire. For case 

(a) the transmitting electron having the same spin as the one in the wire must occupy the 
next level and therefore the energy is increased by eo (by the level separation) . For the case 

(b) when the spins are opposite, the transmitting electron can occupy the same level with 
the electron in the wire. As a consequence, the energy of the system will increase only by 
cqT] (due to the repulsive interaction in the wire). This is the reason for the asymmetric 
resonant condition S{p^^ = = 0) 7^ S{p^^ = 0,pj^ = 1) given in Eq. (15f) and the 
single particle asymmetry in energy, e(p-|- — 1/2 + {r]/2)pi) [see Eq. (12b)]. This formula 
is asymmetric under the transformation — > p^ and p^ — > pi for p^ — 1. As a result the 
current is dominated by electrons which have opposite spins to the one in the wire. It might 
appear from this result that the transmission current is polarized. Since the electron in the 
short wire is in contact with a "thermal bath" the electron has equal probability to be in 
a state \pi — 0,pi — 1) or \p-^ — l,pi — 0) . Due to this fact we conclude that the current 
is not polarized since for each polarized incoming electron only one of the two wire states 
\pi — 0,Pi — 1) or \pi — l,p^ — 0) gives rise to a transmission current. We conclude that 
the conductance for each polarization is half the value of the noninteracting conductance 

Gi e^/2h. Consequently we find that G = G-f + ~ e^/h. 
We remark that Eq. (15g) gives perfect conductance when both e\ and (ei + (77/2 — 1))^ 
obey the condition e^/F^ < 1, (ei + (77/2 — 1))^/F^ < 1. Under this condition we obtain 
from Eq. (15g) G^ ~ e^/^, and consequently G = G| + Gj ~ 2e^ /h. The result obtained for 
this case is in a region where our calculation in Eqs. (12b)-(14b) might not be valid since 
F^ is large. For a short wire 'g^ is large and Eq. (12e) in conjunction with the subsequent 
discussion can invalidate the condition ~ 0. This case is investigated explicitly in the 
next section. 
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V. CONDUCTANCE IN THE LOW TEMPERATURE REGIME 



Equation (10) represents our effective impurity model at lengtfi scale I > la- At this 
length scale the Hamiltonian impurity model himp [see Eq. (Ha)] contains an effective 
Hubbard interaction, ^s(^(i)^/^T^/^J.- finite we can have a situation where the tun- 
neling matrix element A in the coupling Hamiltonian Ht [see Eq. (11c)] obeys A < gs{ld)- 
[A cx cis ~ {d/a)^^'^^''\l/d)]. Here we investigate this case, namely X/cjs < 1- This 

situation is easily achieved in a short wire where gs{ld) is close to the bare Hubbard in- 
teraction [/, which obeys t/U <^ 1. (For a long wire Qs decreases, first exponentially and 
then logarithmically, contrary to A which decreases like dT^/'^. Therefore, for a long wire 
the condition X/qs < 1 is not achieved). In the limit X/gs < 1 we have a strong coupling 
problem. 

The spin degrees of freedom of the impurity interact with the conduction electrons in the 
leads. A virtual charge fluctuation in which an electron migrates off or onto the impurity 
gives rise to a spin-exchange between the impurity and the electrons in the leads. As a 
result an antiferromagnetic interaction between the impurity and the electron in the leads is 
induced [9,20]. Following Eq. (10), we consider only the impurity model plus transmission 
plus the chemical potential part given in Eqs. (9d) and (9h): 

h = himp + hT + h^^^\p^+p^ + rip^p^] - Ea[iPl^{0)7p,,^{0) + ^^^,^(0)4,^(0)] 

+9sik)V^^V^VlV^ + t2X J2 cos(7rp,)[V^t^,(0)V; - V}^o,Ml W 

^=1,1 

where = 0, 1, cr =t, j and < with E^ > 0. 

Using the projection method we project out the double occupancy for Vj, Va and ipl^^{0), 
4,a(0): 

y; = |0)(a| + e,,_,|-a)(T,i|, 

V; = |(7)(0|+e,,_,|T,i)(-(j|, (17a) 

and 

^i. = |0)(<7|+e.,-.|-a)(t,i|, 
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V'o,a=k)(0|+e,,_,|T,i)(-(j|, (176) 

where e-f^j^ — 1 — — ej^^-f. 

The full Hilbert space consists of ly <8) l^(o) = 1 

iy = |0)(0| + |T)(T| + U)ai + |T,i>a,T|, 
i^(o) = |o)(o| + |T)(T| + U)ai + |T,i)a,T|. (17c) 

We project out the double occupancy 

Q^Qv® g^(o) = (I t, T I) ® (I t, i)(i, t I)- (17c?) 
As a result we obtain, in agreement with Refs. [9,20], 

h^ff = -PhQ{QhQ - E)-^QhP = Je//s(0) • S, (18) 

where P = 1 — Q and 

^-(0) = < JO) (I) ^4,/.(o), 'S^=v;t(|) 



Jeff = 2(2A)^ 



1 

+ — 



(19) 



_9s{ld) — Ed Ed_ 

For a short wire such that X/gs{ld) < 1 we can neglect the induced terms h\\ and /ij,. 
(For a long wire this condition is not possible to satisfy.) Therefore, for this case we rely 
on the result given in Eqs. (14e), (15b) and (15e) for single particle transmission. For the 
remaining discussion we consider the strong coupling limit. 

In the strong coupling limit we consider the Hamiltonian obtained in Eq. (18), 

Hk = H^^"^' + JeffS-s(0). (20a) 
The RG equation for the Hamiltonian in Eq. (20a) are [9,20] 

= i"^''^^^' = ^JleffJ±,eff- (206) 
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For J||,e// = J±,eff — Jeff we find [9,20] that Jeff flows to the strong couphng value Jeff{l) — 



i-2j^ffi ■ '^^^^ allows us to introduce the Kondo temperature Tk, 



Tk = — exp 
Kb 



(20c) 



At temperature T < Tk the symmetric electron "orbit" is screened out by the impurity. 
The symmetric screened electron state is given by Cs,a' 

Cs,a = ■^[cLA-d/'^ + CR,a{d/2 + e)] = 2i(sin A;F£)4,a(0). (20o?) 

Here c^^o- and cr^o- ^-re the electrons in the leads. Since the symmetric orbital '0o,(t(O) is 
screened away we project it out into the asymmetric orbital 

Ca,a = -^[cLA-d/^ ' e) ' CR,a(c?/2 + s)] = 2i(sin A;F£)4,a(0). (20e) 

Using the projection operator P = \ipe,a) {^e,a\ we replace and Lo- by PRaP = -^'4'e,cr 
and PLfjP — -^i^ca- We substitute these projections into the boundary term Hbc given in 
Eq. (Ig) and if'««<^« given by Eq. (Ih). Thus we obtain 



He^ yZ ^'"F da;[^|^(x)(-i9a;)^e,a(a^) - (^(2;)(2sinA;i.a)^V'L(2^)^e,a(a^)]- (21) 

Here Hg represents the effective Kondo screened Hamiltonian expressed in terms of the 
asymmetric orbitals ipe,a- For kpa — 7r/2 Eq. (21) represents the single impurity resonant 
level model. This model has perfect transmission, namely G-^ = = e'^/h, and G = 2e^/h. 

In summary, we have obtained perfect transmission, G — 2e^/ h, in the region T < Tk and 
for a short wire such that \/gs{ld) < 1- This is in contrast to a long wire and intermediate 
temperatures, where G ~ e^/h. At low temperatures and a long wire the conductance is 
affected by the two particle transmission, see Eq. (9g). 



VI. DISCUSSION 

To compare our theory with the existing experiments [1,7,8] we used Eq. (15f) to compute 
the conductance. This equation was obtained under the assumption that p^. ~ 0. This 
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assumption is consistent with the assumption of weak transmission; therefore in the Kondo 
regime our formula might not be valid. However, we find that Eq. (15f) works rather well 
and explains both of the regimes of a long wire with weak transmission and a short wire with 
large transmission. Using Eq. (15f) we observe that the shape of the conductance curve is 
mainly dependent on the tunneling matrix element and on the gate voltage Uq = eVa/eQ. We 
plot the conductance G/ {2e^/h), G = G^ + G^ as a function of Uq, for different = (2A)^/eo 
tunneling matrix elements and different P — eo/ZcsT = j'u'feyj' temperatures. 

The conductance is less sensitive to the rest of the parameters 77/2, Qg, electron-electron 
interaction and single particle energy. We use ~ r)/2 ~ 0.75, ei = —0.2 and ic = 1.2 [see 
Eq. (15f)]. In Fig. 1 we plot the calculated conductance as a function of gate voltage Ug for 
fixed transmission F^ = 0.1 (which corresponds to weak transmission) and temperatures (3 — 
1, 10 and 50. At /3 = 1 we observe that the conductance takes the value of G/ {2e^/h) ~ 0.8. 
Decreasing the temperature to /9 = 10 we observe that the conductance has become 1 and 
in addition we observe a shoulder in the conductance around 0.7. At /3 — 50, the shoulder 
develops into a separate conductance peak. Note that this figure very much resembles the 
experimentally measured conductance [1,7,8]. 

To understand the situation for large tunneling matrix elements (i.e. a short wire and in 
the Kondo regime) we chose F^ = 5 and /3 = 100. In Fig. 2 we plot the calculated conduc- 
tance as a function of gate voltage for these parameters. We observe that the conductance 
is G/2e'^/h = 1 as is expected from the Kondo solution given in Sec. V. However, we note 
that in this regime the approximation used in Eq. (15f) is not valid and the Kondo solution 
given in Sec. V should be used. 

VII. CONCLUSION 

We have investigated the case in which a two dimensional electronic waveguide with a 
varying width D{x) in the y direction and e-e interactions can be projected to a one channel 
problem. In this case the reservoirs are replaced by a one dimensional Fermi liquid (leads) 
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coupled to a short wire Luttinger liquid controlled by the gate voltage /z"'*''^. Using a com- 
bined method of Renormalization Group and zero mode bosonization, we have constructed 
the effective spectrum of a short interacting wire. This spectrum consists of combined charge- 
spin density waves and is dominated by the zero mode fermionic spectrum. At temperature 
T < T"'"'^ the fermionic spectrum is equivalent to a multilevel spectrum for which we find the 
conductance formula given in Eq. (15b). Transmission is dominated by the spin S = ±1/2 
doublet giving rise to a conductance G ~ e^/h. This result is relevant for experiments in 
which the wire is short, T < T"'*''^ ~ hvp/kBd. On the other hand, at lower temperatures 
when the effective tunneling matrix A is smaller than the effective electron-electron inter- 
action Qsild), ^/gs{ld) < 1) we find that the short interacting wire is mapped into a Kondo 
problem which gives rise, for T < T^ondo < T"'*''^, to a conductance G — /h. By tuning 
the length of the wire and the temperature we obtain /h < G < 2e^/h in agreement with 
experiments [7]. 
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IX. APPENDIX A 

Here, using the RG equation we derive the coupling constants for the model in Eq. (9a). 
(We use the scaling results for the short wire.) We treat, within a perturbative RG, the 
transmission term Ht- We separate 4>a{x) = (p^i^) + ^(pai^), 4'o,a{^) = '4^o,a{^) + ^'4^o,a{x), 
i^e,(j{^) — '4^e,ai^) ''<' ^'4^e,a{x) ^ud integrate 5(f)a{x), Sipo^ai^), s-nd 5ipe^a-{x) reducing the cutoff 
from A to A' = A — c?A = Ae~K As a result of the integration we replace Ht with HT,eff = 

+ dH^^ + dH^\ where represents the transmission term at the reduced cutoff and 
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fields and dH^^ represents the first order correction in A: 

^ R=±{dl1-e) 

-e-^'=-^xl<(x)V^<,(a;)]. (^1) 

In Eq. (Al) we have used the bosonic representation of the fermions in the leads, ipo,a{x) — 
__i^QiV^Oo,a{x) g^j^^j ^Yie bosonic part (I/VStto) exp{i\/A7i(pa-{x)) of the short wire, see Eq. 
(2c). We obtain 

E «.) - y (^2) 

o-=T,J. 

^ I oTT 

where and Kg are scale dependent (see [8]). Using Eqs. (A2) and (A3) we find 

d\ + + 



dl 



A. (^4) 



To a good approximation Kc + Kg — 2. This gives dX/dl X with the scaling law X{1) — 
Xe~K Therefore X{ld) — X{d/a)~^. Next, we rescale the fields such that '4>o,a{x) is replaced 



by ilJo,a{x), il)o,a{x) = J d/ ailJo,a{x) , where ipo,a{x) depends on the new cutoff A = \/d. Since 



at the scale / > the bosonic field has disappeared, we replace = l/y2TmV^ by V„. As 
a result the new coupling constant becomes, 

1/2 



(7\ i/^ -| , 

- ^= = AA^/^ A=^sinA;Fa, k = l/d. 
a J V27ra V'^ 



{Ah) 



The induced RG terms are obtained from the second order term in A, dH^^ — dH^^p_^ + 

(2) (21 (2) 

dHrpp^p, where dHj.p_i^ represents the particle-hole and dHrp'^_^ the two-particle induced 
interactions: 

dH!^^ = UiXf Y: E [dt[ dTTM<{Q, t+T)e^'-\<{R, t+T)e-^'^-^' xlf{R' , t) 

<T,a'='[,l R,R'=±(d/2-e) 

V'„V(0, t) + e-^-«xt<(i?, ^ + t)V'o1.(0, t + t)<< (0, t)e^'^-^' x<,{E! , t)] [47r(5e„,,(i + t)5V W) 

+AT:{5(t)a{R^ + T)5(t)a'{R!^))] 
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+[C(0' t + r)e''^\^{R, t + r)e''-'''x<,iR', t + r)<<,(0, t) 

+e-''-\i<{R, t + t)^<M t + r)^f,A^, t + T)e-^'-^'xif{R', t)][47r(5^„,,(t + r)5V(^)) 

+A7r{5Mt + r)6(l)a'm}. (AG) 

Using the short wire RG equation we compute the expectation values of {6(1)^6 (pa')- Conse- 
quently, we obtain the RG equations for and /|| which appear in Eqs. (9c)- (9g): 

dl VpJVylj 2 z 



^4 = + + !!£[(MM))l±(MM)):i(l + eosA(Ofl) 

dl " vfA{1) 2 ^ 2 ^ ^ ^ ^ 

_(M)Z^,M±(M)Z^(l_eosA(i)fl)l}. (^8) 



The RG equation for the transversal exchange coupling constant J± is given by 

-(1 -cosA(Oi?)]}. (A9) 



at i'i;'A(<) 2 2 



2 

In Eqs. (A7), (A8) and (A9) A obeys the scaling equation (A4), A(/) = Ae"' and Kc{l), 
Ks{l), Vc{l) and Vs{l) are determined by the sine-Gordon scaling (see [14]). 

The values of J\\{ld)-, J±{ld) and /||(/d) are obtained by solving Eqs. (A7)-(A9) with the 
initial conditions /||(0) = J||(0) = J±{0) = 0. The induced coupling constants are given 
by I\\{ld), J\\{id) and J±{ld)- These couplings are obtained from I\\{ld), J\\{ld) and J±{ld)- 
= /||(Zd)(l/27ra)(d/a), J||(Zd) = J||(/<i)(l/27ra)(d/a), and Mk) = Mk)il/27ra)id/a). 
The values /||(/d), J\\ild) and J±{ld) depend on the sine-Gordon scaling of the short wire. 

Next we consider the two cases: 

(a) {G — AkF)d < 27r. As a result the backward term controls the scahng, Kc/vc ^ 0, 
Ks/vsr^ {1/2vf){1-U/71vf). 
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Using 



we find 



with 



U/a I c dl 

' 



c \ I d\ 



A;27r^l37 \vfJ \a J 



(3-2Kc) 



1 - U/ttIv^ 



l-(l/137)(c/t;^)(l/7rA;) 

\ 1 + U/t:Iv+ + {l/m){c/vF)(l/T:k)' 

Using gc{ld) and Kc we obtain the exchange couphng constants. 



U 



27TVF 



[1 - (cos Ad)) 



hil'i) = P^^og(-] [1 + (cos Ad)], 



(AW) 

{All) 



and j±{ld) is given by 



J±{ld) 



ZTTVp \a 



1 - 



U 

TTVp 



[1 - (cos Ad)) 



(A12) 



In Eqs. (A10)-(A12) A(0) is given by A(0) = 2V2tsm{kFa). U is the original Hubbard 
interaction and (cos AR) is given by 



1 r'-d If 
(cosAi?) = — / cos(Ae ^)dl = — 

Id Jo Id Jh 



1 [R/o- cosx 
ax- 



1 f^/"- . cosx 
ax . 



X 



{A13) 



Id JR/d X log(d/a) Ji 

(b) For the generic case {G — AkF)d > n we can neglect the Umklapp term. We use the 
relations KcVc — KgVs = vp — ^vp] 

Kc_Kl_ 1 (1-U/t^v+\ 



vp 2vf \1 + U/7rv+ J 



and 



Vs{l) 2vf 

As a result we obtain 



2vf 2 Jo 



2vf 



O 



2vf 



, 1 2C/ 2t/ 
~ I 1 h 1 H 



Vc ~'~ Vs 2vf \ 7ri'+ 



nv- ) vf 
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rs^ 



and 



Vc Vs XTTVp/ Vf 

Substituting these values into the scahng equations (A7)-(A9) we find: 



X. APPENDIX B 



Here, we consider the boundary term in the leads with the overlapping region £ ~ a. 

_[(g-ifc^a^t (a) _ e'''^''Rl{-a)){e'''^'R^{e) - e-''''''R^{-e))] 

= to E sin(A;Fa) sm{kFs)[{Ll{-a) - Ll{a)){L^{-e) - L„{e)) 

+ (4 (a) _ i?t _ R^^-e))] + i/.c. (51) 

The leads Hamiltonian in Eq. (le) can be written in terms of even and odd chiral 
fermions: 



The boundary term Hbc becomes, 



Hbc = |(sin(A;^a)) sin(A;^£) ^ {[(V'l.(-a) - - (V'J,.(a) - <.(-a))] 

[(V'e,.(-£) - i^eAe)) - i^oA-e) - ^l^oAe))] + - + - 

[(V'e,a(£) - V'e,a(-£)) + (V'o,.(£) - + He. (S4) 
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The representation (li) is used in the Kondo regime where the symmetric state ipo,a{x) 
is screened out and Hbc is replaced by the antisymmetric state ipe^^^x). In the Kondo hmit 
Hbc is replaced by a "mass" term [the designation symmetric or antisymmetric is given 
according to the original fermions CR^ai^) and cl,(j{x)]-i 

Hbc - -2toi2sm{kFa)y ^ V'U(0)V'e,a(0). (S5) 
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FIGURES 

FIG. 1. Conductance (G) as a function of gate voltage (t/c) for fixed weak transmission 
(r^ = 0.1) for three different temperatures: P = 1 (dotted line), /3 = 10 (dashed line) and /3 = 50 
(solid line) . Note that decreasing the temperature to /3 = 10 leads to the conductance reaching the 
value of 1.0 and in addition there is a shoulder around 0.7. 

FIG. 2. Conductance as a function of gate voltage for large transmission (F^ = 5.0), i.e., for 

a short wire in the Kondo regime at very low temperature (/3 = 100). The conductance is about 
1.0 X 2e^ /%, as expected from the Kondo solution. 
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